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[BAHL DECOPIN GI APPARATUS AND METHOD OF STORINC^ 
REFERENCE VECTOR OF STATE MF.TRTr 



FIELD OF THE INVENTION 



[0001] 



The present invention relates to decoding and in particular to methods 



for implementing Bahl decoding. 

BACKGROUND OF THE INVENTION 
[0002] Messages transmitted over a wireless medium are encoded in order to reduce 
the chances that noise will corrupt the message. The encoding translates an original 
block of an original message which comprises a plurality of bits which are either '0' or 
•1' to an encoded block which is formed of a plurality of symbols which are either '0' 
or T, The encoding of the block usually includes adding redundancy (i.e., making the 
encoded block longer than the original block), introducing symbol dependency (i.e., 
making the value of each symbol in the encoded block a function of a plurality of bits 
in the original block) and interleaving the block (i.e., mixing the order of the symbols 
in the encoded block). In modem codes the size of the original bloclc is usually very 
large, of the order of hundreds, thousands or tens of thousands of bits. Thus, tihe 
encoding turns a block of binary bit$ into a block of binary symbols. It is noted that 
there is not necessarily a direct correlation betweei the value of a single bit and the 
value of any one of the encoded symbols. However, the number of bits in the original 
block and the number of symbol groups in the encoded block are generally the same. 
Generally in convolutional codes, bits and symbol groups in coiresponding locations 
within the block are related in that they form pairs of input and output of encoders aad 
decoders, as described hereinbelow. 

[0003] In convolutional codes, an encoder which has a fixed number of possible 
states is used to encode the original block. Each bit of the original block which enters 
the encoder changes a current state of the encode to a next state according to the 
value of the entering bit. It is noted that for convolutional codes, from each state the 
encoder can move only to two other states based on whether the input is '0' or '1'. For 
each entering bit the encoder generates a group of one or more symbols forming the 
transmitted message. In other trellis codes, larger inputs are received by the encoder 
in each iteration, and the number of possible state changes are larger accordingly. In 
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decoding, the decoder assigns each state a. state metric (SM), which represents the 
conditional probability that the encoder was in that state. 

[0004] During transmission of the message, noise is added to the transmitted encoded 
message, and a decoder is used to extract the original message from the received 
message, which includes the noise. Normally, each symbol in the received message is 
represented in the receiver, before decoding, by a word which represents the chance 
that the symbol is a logical 1'. A high positive value of the word means that the 
symbol is a logical '1' with high confidence, while a negative value with a high 
absolute value means the symbol is a logical '0' with high confidence. A word with a 
zero value means that the symbol is a '0' or 'V with equal probability. This 
representation is referred to herein as soft data. 

[0005] Some decoding methods include one or more iterations of a step of an A- 
Posteriori probability estimation in which the soft data representation of each symbol 
is adjusted and/or represented by a hard value (i.e., either '0' or 'l*) based on the values 
of other symbols in the block of the adjusted symbol. 

[0006] Many methods for perfonning the adjustment step are based on the Bahl 
algorithm which is described^ for example, in Steven S. Pietrobon, "Implementation 
and Performance of a Turbo/MAP Decoder", Int. J. Satellite Commun., vol. 16, pp. 
23-46, Jan, -Feb. 1998, the disclosure of which is incorporated herein by reference. In 
the Bahl algorithm, for each received symbol group, two vectors of state metrics are 
prepared: a forward vector which represents the conditional probability of the encoder 
being in each of the possible states immediately before the symbol group was 
generated and a backward vector which represents the conditional probabiUty of the 
encoder being in each of the possible states immediately after the symbol group was 
generated. The forward vector is calculated based on the received symbols which 
precede the adjusted symbol group in the block and the backward veaor is calculated 
based on the received symbols which follow the adjusted symbol group in the block. 
The forward and backward vectors are recursively dependent on the previous state 
metric vector in the respective directions from the end points of the block. 
[0007] The forward and backward vertors of each symbol group are suppUed 
concurrently to a log likelihood ratio (LLR) unit which provides an adjusted output 
value of the symbol group based on the forward and backward vectors and the 
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Original value of the symbol group. The output value is chosen according to the most 
likely transitions between the state metrics of the forward and backward vectors. 
[0008] After the vectors are supplied to the LLR unit, the vectors may be discarded 
unless they are needed for calculating vectors of other symbol groups. Therefore, it is 
preferred, in order to save storage, that the vectors be calculated immediately before 
they are used. However, due to the recursive nature of the vertors which are 
calculated from opposite sides of the block, such immediate calculation, usage and 
discarding is performed m the prior art only in one direction, e.g.. in the forward 
direction. 

(0009] A conventional prior art method for calculating the output values of all the' 
symbols in a block comprises two passes over the block. In a first pass over the block, 
performed before the LLR unit begins the adjusting, the backward recursion vectors 
for all the symbols in the block are calculated and stored. In a second pass over the 
block, beginning from the opposite end of the block, the forward recursion vectors are 
calculated, the respective backward recursion vectors are retrieved from the storage 
and the forward and backward vectors are supplied to the LLR unit. The supplied 
vectors are immediately used by the LLR unit and then they are discarded. 
[0010] The storage requirements in this method, required in order to store the 
backward recursion vectors, are quite large. In order to reduce the required storage 
space, it has been suggested to perform the symbol adjustments in sub-blocks in 
which each symbol is influenced only by the symbols in its sub-block, rather than by 
all- the symbols in the block. Such method is described, for example, in PCT 
publication WO98/20617 the disclosure of which is incorporated herein by reference. 
The WO98/20617 publication shows that by beginning the recursion calculations 
from a point sufficiently before the beginning of the sub-blodc. the recursion values 
corresponding to the symbols in the sub-block are sufificiently valid although the 
recursion did not start from the b^inning of the block. Even so, the storage 
requirements for a 16-state code, using 10-bit state metrics for blocks of 40 symbols, 
is 6400 bits. When implementing a hardware decoder such an amount of storage is 
very expensive. 

[001 1] Fig. lA is an exemplary trellis map 10 of the convolutional code 7.5 which is 
known in the art. The encoder of the 7.5 code may be in any of four states s(n) 12 
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identified by a combination of two bits. From each of the states s(n), the encoder can 
move to one of two other states s(n+l) based dn an input bit d(n), as indicated by 
branches 14. In moving from state s(n) to state s(n+l) 16, the encoder provides a 
symbol group formed of a pair of bits co(n) and ci(n). The symbol group co(n) and 
ci(n) is transmitted to a receive which uses a decoder to determine the original bits 
d(n). The decoder assigns each received symbol group with a vector of branch metric 
(BM) values 6n(co,Ci) which designate the conditional probability that the received 
symbol group originates from the transmitted symbol group designating cq and ci. In 
addition, the decoder calculates forward (an) and baclcward (pn) state metric vectors 
which have elements for each of the states s(n) which designate the conditional 
probability that the decoder was in the state s(n) at time n. Generally, forward vectors 
an are recursively dependent on preceding vectors (on+i = f(a„ , Sq)), and backward 
vectors Pn are recursively dependent on following vectors (P^ = f(3n+l . 5n))- 
[0012] As can be seen in Fig. 1 A, Each pair of states s(n) is connected to a pair of 
states s(n+l) although not necessarily the same states. The pairs of states with the 
connecting branches form a butterfly trellis which is depicted in Fig. IB which is a 
butterfly trellis diagram for convolutional codes. This feature of pairs of states is 
generally true in trellises of convolutional codes. In Fig. lA states 00 and 01 (which 
are equivalent to states X and Y in fig. IB) at time n are connected to states 00 and 10 
(which are equivalent to states Z and W in fig. IB) at time n+1. Similarly, states 10 
and 1 1 at time n (which are equivalent to states X and Y in fig. IB) are connected to 
states 01 and 11 at time n+1 (equivalent to states Z and W in fig. IB). 

SUMMARY OF THE INVENTION 
[0013] An aspect of some preferred embodiments of the invention relates to reducing 
the amount of memory required for decoding convolutional codes. 
[0014] An aspect of some preferred embodiments of the present invention relates to a 
method of providing pairs of forward and baclcward vectors concurrently to an LLR 
unit for A-Posteriori probability estimation of a decoded block. During a first stage, 
before the beginning of the operation of the LLR unit, substantially all of the 
backward vectors are calculated but only one or more reference vectore, from which 
the rest of the vectors may be easily re-calculated, are stored. During a second stage, 
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the backward vectors which were not stored are calculated from the stored reference 
vectors and are supplied to the LLR unit along with the forward vectors. 
[0015] In a preferred embodiment of the present invention, during the first stage, 
reference vectors are stored at predetermined intervals. In the second stage the 
reference vectors are retrieved in course from the storage and are used to re-calculate 
a group of non-stored vectors between a current and a previous reference vector. 
[Preferably,] A short term storage area is used to store the re-calculated vectors. Each 
group of non-stored vectors are preferably stored in place of the previously re- 
calculated vectors, as the previously re-calculated vectors were used ateady and are 
not needed any more. 

[0016] An aspect of some preferred embodiments of the present invention relates to 
using a reverse fiinction in recalculating at least some of the vectors. In some 
preferred embodiments of the invention in which the backward vectors are calculated 
in the first stage using a reversible jEunction, the backward vectors are calculated in the 
second stage using the reverse fonction of the function of the first stage. Thus, both 
the backward and forward vectors are calculated in the same direction of progression 
over the block during the second stage. In a preferred embodiment, the vectors 
calculated in the first stage are not stored. Alternatively, a plurality of backward 
vectors from the first stage are stored, in order to limit possible round ofif errors due to 
the use of the reverse fiinction and/or when singular functions are encountered. 
Further alternatively or additionally, the reverse function is used to re-calculate only 
some of the non-stored vectors and the rest of the recalculated vectors are calculated 
using the non-reversed function. 

[0017] [There is therefore] Therefore there i<! provided in accordance with a preferred 
embodiment of the present invention, a method of calculating recursive state metric 
vectors of a block of symbols, the state metric vectors being supplied to an output 
calculating unit, including pre-calculating the state metric vectors from a first side of 
the bloclc, and calculating for a second time at least some of the state metric vectors 
from the first side of the block. 

[0018] Preferably, pre-calculating the state metric vectors includes pre-calculating the 
state metric vectors before supplying substantially any of the vectors to the output 
calculating unit. 
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[0019] Preferably, the method includes storing fewer than all the pre-calculated 
vectors. 

(0020] Preferably, calculating for a second time at least some of the state metric 
vectors includes calculating those vectors which were not stored. 
[0021] Preferably, storing fewer than all the pre-calculated vectors includes storing a 
number of vectors which is about the square root of the number of pre-calculated 
vectors. Preferably, storing fewer than all the pre-calculated vectors includes storing 
vectors selected responsive to the locations of singular functions used in the pre- 
calculating of the vectors. Alternatively or additionally, storing fewer than all the pre- 
calculated vectors includes storing vectors selected in predetermined intervals. 
Preferably, storing vectors seleaed in predetermined intervals includes storing vectors 
with equal intervals between them. Alternatively, storing vectors selected in 
predetermined intervals includes storing vectors with intervals of decreasing size 
between them, 

[0022] Preferably, calculating for a second time the at least some of the state metric 
vectors includes calculating the vectors based on the stored vectors. 
[0023] Preferably, calculating for a second time the at least some of the state metric 
vectors includes calculating at least some of the at least some of the vectors using a 
reverse function of a function used in the pre-calculating. Preferably, substantially 
none of the pre-calculated vectors are stored during the pre-calculation. Alternatively, 
substantially only pre-calculated vectors from which a new pre-calculated vector is 
calculated using a non-reversible flinction, are stored during the pre-calculation. 
[0024] Preferably, calculating for tiie second time includes calculating at least some 
of the vectors using the funaion which was used in the pre-calculation. Alternatively 
or additionally, the pre-calculatmg is perfonned using a function which is an 
approximation of an original function and calculating for a second time tiie at least 
some of tiie state metric vectors includes calculating at least some of the at least some 
of the vectors using a reverse function of the original fiinction. Further alternatively or 
additionally, the pre-calculating is performed using a function which is an 
approximation of an original function when tiie original function is non-reversible. 
[0025] Preferably, pre-calculating tiie state metric vectors from a first side of tiie 
block includes calculating forward state metric vectors. Alternatively, pre-calculating 
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the state metric vectors from a first side of the block includes calculating badcward 
state metric vectors. Preferably, the method includes calculating state metric vectors 
from a second side of the block after the pre-calculating. 

[0026] Preferably, pre-calculating the state metric vectors includes pre-calculating a 
number of vectors substantially of the size of an encoding block. Alternatively or 
additionally, pre-calculating the state metric vectors includes pre-calculating a number 
of vectors substantially smaller than the size of an encoding block. 
[0027] Tliere is further provided in accordance with a preferred embodiment of the 
present invention, a method of calculating recursive state metric vectors from opposite 
sides of a block of symbols, the state metric vectors being supplied to an output 
calculating unit, including pre-calculating a plurality of state metric vectors from a 
first side of the block, storing fewer than all the pre-calculated state metric vectors, 
and calculating state metric vectors from a second side of the block after pie- 
calculating the plurality of state metric vectors from the first side of the block. 
Preferably, the method includes calculating at least one of the state metric vectors 
from the first side of the block for a second time after calculating at least one of the 
state metric vectors from the second side of the block. Preferably, calculating at least 
one of the state metric vectors from the first side of the block includes calculating the 
state metric veaors from the stored vectors. Preferably, storing fewer than all the pre- 
calculated vectors includes storing vectors calculated in predetermined intervals. 
Preferably, calculating the at least one of the state metric vectors includes calculating 
based on the closest stored vector. 

[0028] Preferably, the method includes presenting an output calculating unit with 
pairs of corresponding state metric vectors from opposite sides of the block. 
Preferably, the block of symbols is divided into a pluraUty of segments defined by the 
stored vectors and the method includes calculating for substantiaUy all the segments at 
least one of the vectors for a second time based on the respective stored vectors of the 
segment. Preferably, calculating for substantially all the segments includes calculating 
for aU the segments except the first segment. Preferably, the method includes storing 
the at least one of the vectors which was calculated for the second time. Preferably, 
calculatbg the at least one of the vectors for the second time includes calculating 
substantially all the non-stored vectors of the segment. Prefefably, the method 
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includes storing substantially all the vectors of the segment which were calculated for 
the second time. Preferably, the method includes calculating at least one of the vectors 
from the first side of the block for at least a third time. 

[0029] There is fiirther provided in accordance with a preferred embodiment of the 
present invention, a method of calculating recursive state metric vectors of a block of 
symbols, the state metric vectors being supplied to an output calculating unit, 
including calculating a plurality of state metric vectors representing a first side of the 
block using a first fiinction, and calculating a plurality of state metric vectors 
representing the first side of the block using a second function different fi-om the first 
function. 

[0030] Preferably, the first function is identical to a third function when the 

third function is reversible and is received by adding a small vahie to the third 
fiinction when the third function is non-reversible. Preferably, the second function is 
substantially a reverse of the first funaion. 

[0031] Alternatively, the first fiinction includes an approximation of a third function 
and the second function includes a substantial reverse of the third function. 
[0032] Preferably, the first function includes a log MAP function. 
[0033] Preferably, calculating the plurality of vectors using the first fimction includes 
calculating substantially all the vectors representing the first side of the block. 
[0034] There is fijrther provided in accordance with a preferred embodiment of the 
present invention, a method of calculating recursive state metric vectors of a block of 
symbols, the state metric vectors being supplied to an output calculating unit, 
including calculating a plurality of state metric vectors representing a first side of the 
block recursively based on immediately preceding vectors in the block, and 
calculating a plurality of state metric vectors representing a second side of the block 
based on immediately preceding vectors in the block, 

[0035] Preferably, calculating the plurality of vectors representing the first and 
second sides of the block includes calculating pairs of vectors representing first and 
second sides of the block, substantially concurrently. 

[0036] Tliere is further provided in accordance with a preferred embodiment of the 
present invention, a state metric calculation unit, including circuitry for recursively 
calculating state metric vectors from a first side of a block including a predetermined 
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number of symbol groups, and a memory for storing state metric vectors which has a 
maximal storage space for state metric vectors substantially smaller than that required 
to store the predetermined number of state metric vectors. 

[0037] Preferably, the.maximal storage space of the memory has room for storing less 
than fifty percent of the predetermined number of state metric vectors. FurthM- 
preferably, the maximal storage space of the memory has room for storing less than 
twenty percent of the predetermined number of state metric vectors, 
[0038] Preferably, the circuitry implements a plurality of different functions for 
caJculating the state metric vectors. Preferably, the circuitry implements a pair of 
functions for calculating the state metric vectors which pair include mutual reverse 
functions. 

[0039] Preferably, the memory includes a long term storage area for storing pre- 
calculated state metric reference vectors in predetermined intervals and a short term 
storage area for storing state metric vectors between two reference vectors. 
Preferably, the long term storage area serves for storing also state metric vectors from 
between two refwence vectors. 

BRIEF DESCRIPTION OF FIGURES 
[0040] The invention will be more clearly understood by reference to the 

foUowing description of preferred embodiments thereof in conjunction with the 
figures, in which same number designations are maintained throughout the figures for 
each element and in which: 

[0041] Fig. lA is an exemplary trellis map of the convolutional code 7.5 which is 
known in the art; 

[0042] Fig. IB is a butterfly trellis diagram for convolutional codes which is known 
in the art; 

[0043] Fig. 2 is a schematic illustration of a decoder, in accordance with a preferred 
embodiment of the present invention; 

[0044] Figs. 3A and 3B are a flow chart of the actions performed in decoding 
received data symbols, in accordance with a preferred embodiment of the present 
invention; 

[0045] Fig. 4 is a schematic illustration of SM vectors of a decoded sub-block, in 
accordance with a preferred embodiment of the present invention; 
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[0046] Fig. 5 is a schematic illustration of SM veaors of a decoded sub-block, in 
accordance with another preferred embodiment of the present invention; and 
[0047] Fig. 6 is a schematic illustration of SM vectors of a large block, in accordance 
with a preferred embodiment of the present inventioa 

DETAILED DESCRIPTION OF THE PREFERjRED EMBODIMENTS 
[0048] Fig. 2 is a schematic illustration of a decoder 20, in accordance with a 
preferred embodiment of the present invention. Decoder 20 preferably receives soft 
data symbols along an input channel 22. A buffer memory 24 accumulates a sufRcient 
number of symbols in order to decode a sub-block. The received symbols are grouped 
into groups of one or more symbols according to the encoding scheme of the received 
data symbols. For each symbol group, an LLR unit 30, calculates an output value 
corresponding to the location of the symbol group in the sub-block. The calculation is 
preferably based on the values of substantially all the symbol groups in the sub-block. 
35 is known in the art. and described, for example, in the above mentioned publication 
of Steven S. Pietrobon. In the following description the accumulated sub-block is 
assumed to comprise LI groups of symbols (Ri,...,Rli). In addition, buffer memory 

24 preferably accumulates L2 additional groups of symbols (Rl1+1 Rl1+L2) 

purpose of which is described hereinbelow. 

[0049] As is known in the art, a branch metric rB TVH calculation unit 25 calculates 
branch metrics (BMs) which represent the groups of symbols from buffer memory 24 
in the calculations. Alternatively, buffer memory 24 accumulates the BMs which were 
calculated by BM calculation unit 25. 

[0050] A forward calculation unit 26 and a backward calculation unit 28 calculate 
forward and backward state metric (SM) vectors, respectively, based on the 
accumulated symbols, as is known in the art. 

[005 1] LLR unit 30 uses the SM vectors from forward and backward calciili^ inn units 
26 and 28, together with the BMs from BM calculation unit 25, to provide the output 
values for each of the groups of symbols in the sub-block, using any method known in 
the art. The output provided by LLR unit 30 for each group of symbols, may comprise 
hard-data and/or soft-data values, may comprise a decoded value or may comprise an 
intermediate value used in additional decoding steps. 

10 
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[0052] In order to calculate its output for each symbol group. IXR unit 30 receives 
the respective BM (5n) of the bit from buffer memory 24 and the SM vectors (oq, 
Pn+l) which correspond to the location of the symbol group within the sub-block, 
from both forward calculation unit 26 (ocn) and baclcward calculation unit 28 (Pn+l)- 
The forward SM vectors (ctn) are calculated by forward calculation unit 26 
recursively from the beginning of the sub-block, based on a starting SM vector (ai) 
and the BMs (5i,....5n.i). The starting SM vector (aj) is preferably determined either 
based on the SM vector from a previous sub-block or based on a predetermined value 
from which the decoding begins, e.g.. a vector of zeros. Likewise, the backward SM 
veaors (Pn) are calculated by backward calculation unit 28 recursively from the end 
of the sub-block based on a last backward SM vector (3li+i) of the sub-block and 
the BMs (5n+i,...,5Li). The ending backward SM vector (Pli+i) is of a 
predetermined value or is calculated from the additional groups of symbols 
(Rli+1>-.Rli+L2) using any method known in the art, for ©cample, that described 
in the above mentioned PCI publication WO98/20617. The number of additional 
groups of symbols L2 depends, as described in the above publication, on the size of 
the encoder and on the desired accuracy of the calculations. 

[0053] Backward calculation unit 28 preferably calculates the SM vectors in two 
stages. In a first preparation stage, performed before the beginning of the operation of 
LLR unit 30, backward calculation unit 28 calculates the backward SM vectors 
recursively and stores a plurality of them as reference vectors to be used in the second 
stage. The reference vectors are preferably stored in a long terni storage 38 of m 
internal memory 34. In the second stage, a decoding stage, backward calculation unit 
28 preferably operates in parallel with LLR unit 30 and forward calculation unit 26 
and calculates the backward SM vector (3n) presently needed by the LLR unit, based 
on the closest stored reference vertor of a higher index (the backward vectors (Pn) are 
calculated from the highest index to the lowest index recursively, using methods 
known in the art). In a preferred embodiment of the invention, internal memory 34 
includes a short term storage area 36 in which the backward SM vectors between two 
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reference vectors are stored between their calculation and use during the decoding 
stage. This is a recalculation of vahies previously calculated and thrown away. 
[0054] A controller 32 controls the timing and operation of the units of decoder 20. in 
accprdance with the following description. 

[0055] Figs. 3 [A and 3B are] is a flow chart of the actions performed in decoding 
received data symbols by decoder 20, in accordance with a preferred embodiment of 
the present invention. Buffer mgmorv [Memory] 24 accumulates (block 50) symbols 
belonging to a sub-block of LI symbol groups. In addition, buffer memory 24 
preferably accumulates (bloc^c 52) L2 additional symbol groups in order to calculate 
the last backward SM vector of the sub-block. BM calculation unit 25 

preferably creates (block 54) respective representative BMs from the accumulated 
symbols, as is known in the art. 

[0056] Reference is also made to Fig. 4 which is a schematic illustration of the 
backward SM vectors (Pn) of the decoded sub-block, in accordance with a prefeired 
embodiment of the present invention. In the preparation stage, backward calculation 
unit 28 calculates (blodc 57) backward SM vectors (Pn), using the representative BMs 
of the accumulated symbols of the sub block. Backward calculation unit 28 first 
determines (IslQck 56) the last backward SM vector (3li+i) 100 of the sub block. 
Preferably, baclcward c^cuiatio^ unit 28 chooses an arbitrary SM vector as the value 
of PLl-f-L2, for example a zero vector. Using the arbitrary vector and the L2 
additional symbol groups, backward calculatini^ unit 28 calculates recursively the SM 
vectors PL2+L1 to Plj+i. The calculation is performed as is known in the art. for 
example as described in any of the above referenced application, With the calculation 
of each of the vectors, the previous vector is [preferably] discarded as these vectors do 
not carry reliable data and are not needed for further calculations. Vector Pli+1 on 
the other hand is considered reliable as it is the result of L2 consecutive symbol group 
calculations and is [preferably] stored in long term storage 38. for later use. 
[0057] After baclcward SM vector is calculated, backward calculatig n unit 28 
[preferably] continues to recursively calculate (block 58) the SM vectors 3li. 
and periodically stores (bJedc 62) some of the vectors 102 in long term storage 
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38, for later use. The reference vectors 102 to be stored are [preferably] chosen 
according to a predetennined scheme. In some preferred embodiments of the 
invention, the reference vectors are distributed throughout the sub block in equal 
intervals of SM vectors 104. In a preferred embodiment, illustrated by blocks 60, 61 
and 64, backward calculation unit 28 stores every ith backward SM vector, i.e., 
PLl+l-mi (ra=O...MsLl/i) where i is close to the square root of LI, e.g., between 6 
and 10. In the following description it is assumed for simpUcity that Ll is divisible by 
i, and the stored backward reference vectors 102 are pj+j, 32i+l. 3l1+1 
[0058] After the last backward vector 102 (pi+j) to be stored in long term storage 38 
is calculated, the decoding stage [Q^OD] is performed. The decoding . stage is 
[preferably] formed of a plurality of steps (M steps). In each step LUR unit 30 
provides output for i groups of symbols. Each of the steps, is [preferably] formed of 
two phases. In a first phase (block 66), backward calculation unit 28 [preferably] 
calculates the backward SM vectors (3mi.....P(m.l)i+l) needed for that step, and 
stores them in short term storage 36. Backward calculation unit 28 [preferably] 
retrieves a reference vector (p^i) 102 from long term storage 38. and calculates 
therefrom, recursively. SM vectors (Elmi-l....J(m.l)i+l). In a preferred embodiment, 
also the retrieved reference vector (p^j) is stored in short term storage 36. 
[0059] In the second phase (bb^ 71) of each step, forward calculation unit 26 
calculates (blQsk 68) the forward SM vectors (On) of that step, and forward and 
backward calculation units 26 and 28 provide (block 72) concurrently pairs of forward 
and backward SM vectors (a^, p„+i) to LLR unit 30. Preferably, backward 
calculation unit 28 provides the SM vectors by retrieving them from short term 
[memory] siorage 36. LLR. unit 30 immediately calculates (block 74) its output 
responsive to the provided vectors and the relevant BM. 

[0060] Decoder 20 then [preferably] proceeds in decoding a foUowing sub-block 
(returning to block 50), as is known in the ait. It is noted that some of the operations 
in handling the following sub-block may be performed concurrently with the 
decoding stage 70 of the current sub-block. 
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[0061] The number of reference vectors stored in long term storage 38 is [preferably] 
chosen in a manner which minimizes the required size of internal memory 34. 
[Preferably,] Jhe minimization is performed in a manner which does not increase the 
required computational resources of backward calculation unit 28 above a desired 
threshold. In a prefeued embodiment of the invention, the number of reference 
vectors i stored in long term storage 38, is approximately equal to the square root of 
LI (the size of the decoded sub-block), as using the square root minimizes the size of 
memory 34. Thus, the long term storage 38 is of a size which accommodates j vectors 
and the size of short term storage 36 is, as well, one which accommodates i vectors. 
Thus, the total size of memory 34 is one which accommodates 2i SM vectors. 
[0062] Alternatively, long term stoi^ 38 is larger, in order to reduce the amount of 
calculations required during the first phase fblock 66) of the steps of the decoding 
stage. Further alternatively, the number of vectors 104 between each pair of reference 
vectors 102 is chosen in a manner which simplifies the implementation of backward 
paloulatiQp unit 28. In a preferred embodimem of the invention, the distance between 
consecutive reference vectors 102 is equal to 2^, n being an integer. Thus, the 
hardware implementation of backward calculation unit 28 is simplified. 
[0063] Alternatively to having the same number of SM veaors 104 between the 
stored reference vectors 102. the number of vectors 104 between the stored reference 
vectors 102 differ throughout the sub-block. In a preferred embodiment of the present 
invention, the number of vectors 104 between the stored reference vectors 102 
increases with the indexes n of the symbol groups of the sub-block, i.e., the number of 
vectors 104 decreases with the progression of the backward calculation fi-om the end 
of the sub-block to the beginning. During the first step of the decoding stage, long 
term storage 38 is full with SM vectors required at later steps of the decoding stage 
and therefore short term storage 36 must be used to store all the SM vectors calculated 
during the first phase (Mock 66) of the first step of the decoding stage. At later steps, 
however, the number of spaces in long term storage 38 which are used decreases and 
the vacam spaces are [preferably] used for storing some of the vectors calculated 
during the first phase (bl.oclc 66) of later steps of the decoding stage. In this 
embodiment, the border between long term storage 38 and short term storage 36 
moves with the steps of the decoding stage [(]70[)]. 
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[0064] For example, in a preferred embodiment of the invention, the size Ll of the 
sub-block is 100. and vectors pio, 320. P31, 343. 356. 370. Pss. 3l01 are stored 
in long term storage 38 during the preparation stage ( block 57). In the first phase 
(bledc 66) of the first step of the decoding stage, backward calculation unit "28 
calculates eight vectors, i.e., 32, - , P9. ""^hich are stored in short term storage 36. In 
the first phase of the second step of the decoding stage, backward calculation unit 28 

calculates nine vectors, i.e., 3n 3ig, 319. Eight of the vectors, e.g., 3n, ... , 3i8, 

are stored in short term storage 36, and the ninth vector, e.g., 319, is stored in long 
tenri storage 38 in place of vector pjQ which is not needed anymore. In each 
following step, the number of vectors calculated in the first phase (block 66) increases 
by one, and the additional vector is stored in long term storage 38 in the places of the 
vectors which are not needed anymore. Thus, the number of spaces required for 
storing vectors in memory 34 is equal to 17, relative to 19 required spaces for evenly 
spaced reference vectors. 

[0065] For some decoding methods, such as for the log-map method, a fiinction f 
used in calculating backward SM veaor Scorn (3n+l) is a reversible fiinction, 
such that a reverse function f-1 may be used to calculate (3n+i) from (3n). In some 

preferred embodiments of the present invention, the reverse function f-1 is used to 
calculate the backward SM vectors (3ji) fi-om the reference vectors during the 
decoding stage [G700]- 

[0066] In a preferred embodiment of the present invention, the first backward SM 
vector of the sub-block (32) is calculated during the preparation stage (block 57). 
During the decoding stage, backward calculation unit 28 calculates vector (3n+i) 
from (3n) using f^l, concurrently with the calculation of (a^) from (0^.1) by forward 
calculation unit 26. [Preferably, in] In this embodiment memory 34 comprises only 
long term storage 38, and the first phase (block 66) of the decoding stage [QTOD], 
which is time consuming, is eliminated. 

[0067] In a preferred embodiment of the present invention, backward calculation unit 
28 does not include aji internal memory 34'. Rather, during the preparation stage 
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(hlQcls 57) backward calculation unit 28 calculates the first backward SM vector O2) 
and uses the reverse function f-^ to calculate all the backward SM vectors jSrom 32, 
during the decoding stage [Q 70D]. 

[0068] It is noted, however, that in some cases fl introduces inaccuracies in the 
calculations due to round off problems of reverse functions. Therefore, as described 
above, in some preferred embodiments, reference vectors are stored in long term 
storage 38 during the preparation stage (^[ock 57), and the reverse function f^l is used 
only between the reference vectors. [Preferably, the] Tlie distance between the 
reference vectors is set based on a maximal number of allowed consecutive reverse- 
function iterations R. The maximal number of reverse iterations R is [preferably] 
determined so that the round off error due to the use of the reverse function f^l is 
limited below a desired level. 

[0069] Fig. 5 is a schematic illustration of the SM vectors (Pn) of the decoded sub- 
block, in accordance with a preferred embodiment of the present invention. In some 
preferred embodiments of the present invention, only some of the backward SM 
vectors (Pn) are calculated using the reverse function f^l during the decoding stage 
[(]70[)], The rest of the backward SM vectors (3n) are [preferably] calculated during 
^ the first phase (block 66) of the decoding stage [0700] using the forward fiinction ^ 

as described above. For example, for each group of SM veaors 104, the first R 
I vectors 107, which have lower indexes are calculated using f^l. The rest of the 

! vectors, i.e., vectors 109. are [preferably] calculated from the nearest reference vector 

I 102 using function f, as described above. Stated otherwise, in each step of the 

decoding stage [(]70[)] backward calculation unit 28 retrieves a reference vector 102 
from long term storage 38. Using fimction f, vectors 109 are calculated and stored in 
short term storage 36, during the first phase (blosk 66). Thereafter, the same reference 
vector 102 is used with the reverse funttion f-1 to calculate SM vectors 107. 
[0070] In a preferred embodiment, the distance between the reference vectors is 
determined using any of the considerations described above and ba^ed on the maximal 
number of allowed consecutive reverse iterations R Designating the number of 
vectors which may be stored in short term storage 36 by Si, the number of vectors 
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between two neighboring reference vectors 102 is Sl+R. The number of reference 
vectors is S2=L1/(S1+R) and it is desired to minimize the space required in memory 
34, i.e.. S1+S2. 

[0071] In some preferred embodiments of the present invention, funaion f is a non- 
reversible function fj which is an approximation of a reversible ftinction f2, e.g., due 
to round off or when using the log-max method which is an approximation of the log- 
map method. In one of these embodiments, some or all of the calculations which are 
performed using f; are performed with Mj, the non-reversible function. Calculations 
performed with the reverse function f-l and possibly some of the calculations 
performed using f, are performed with fj (f2"l). 

[0072] FoUoAving is a showing of the reverse function f-1 for exemplary functions f 
Generating the reverse function equations for any other code will be obvious to those 
skilled in the art. 

In the 7.5 code described hereinabove with reference to fig. lA, using the 
linear MAP method, fiinction f is defined by the equation set: 



11 



3? 

<f 1 1 oOO 



and 



Ai 
.A 



10 



ill 



+ 



(1) 



The reverse function is therefore defined by the following equations: 



cll rOO 
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and 
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[0073] In some preferred embodiments of the present invention, when it is determined 
during calculation of a veaor Pq fi-om 3^+1 in preparation stage (block 57) that the 
branch metrics 6n form a non-reversible (singular) function, one or more of the 
branch metrics 6n are slightly changed in order to remove the singularity. For 
example, referring to equations (1), when for a specific n. the matrix in either of 
equations (1) is singular, i.e., a value le for a small e is 

added to the matrix. That is, a value s is added to l^P^ or b^^^ as required. 

Alternatively, a value e is added to Sn^l or b^}^ as required in order to remove the 
singularity. 

17 



Received from < 972 9 9574933 > at SI5I03 8:35:10 AM [Eastern Dayliglit Time] 



5. MAY. 2003 14:39 EITfiN^EfiRLR LATZER+COHEN-ZEDEK . on275 P. 6 



[0074] AJtematively or additionally, when the fiuiction for calculating vector pj, from 
vector 3n+l is singular, [calculating] backward calculation unit 28 stores p^+i such 
that it is not necessary to use the reverse of the singular matrix. [Preferably, memory] 
Memorf 34 comprises a special storage area (not shown) for vectors stored due to 
singularities. Alternatively or additionally, at least some of the reference vectors 102 
are chosen based on the singularities. In a preferred embodiment of tiie present 
invention, bacl<;ward calculation unit 28 stores reference vectors at equal intervals 
unless a singular equation is encountered in which case the. vector on which tiie 
singular fiinction is applied is stored and the equal vector intervals are [preferably] 
counted from that vector. In another preferred embodiment of tiie invention, the 
number of reference vectors 102 stored during tiie preparation stage fblock 57) is 
constant and the intervals between reference vectors 102 when singularities are not 
encountered are adjusted according to tiie number of Singularities encountered during 
tiie preparation stage. For example, preparation stage fblock 57) may begin with 
storing each tentii veaor. If a large number of singularities are encountered the 
interval between reference vectors is enlarged and. for example, every 15th veaor is 
stored. If. however, less than an average number of singularities are encountered, the 
interval between stored vectors may be shortaied. 

[0075] Ahematively to storing entire vectors to which singular ftmctions are 
applied, a part of these vectors 0^+1. including one or more elements, is stored. The 
one or more stored elements [preferably] allow calculation of tiie rest of tiie elements 
Pn+1 frora Pn despite the singularity of the fimction. 

[0076] AJtematively or additionally to storing vectors to which singular fiinctions are 
applied, vectors which cannot be calculated using the reverse fiinction f-1 due to tiie 
singularity, are calculated during tiie first phase ( block 66) of the decoding stage [Q 
70[)] using function f [Preferably, for] For each interval between reference vectors 
102, the number of vectors 109 (Fig. 5) calculated using function f; and the number of 
vectors 107 calculated using fimction f-^, are dynamically adjusted according to 
whether the interval includes singular functions and where tiie singular fimction 
appears within the interval. 
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[0077] For the log MAP method, ^ ^^^^ values which are the logarithms of 
their values in the linear case. Using the definitions in Fig. IB, the equations defining 
f for a pair of states x and y, are: 

= logie^'^-^^^l + ) and = log(e^^^^*l + e"^' +^"+1 ) . 

The revCTse function f-1 is [preferably] defined by the equations: 



yCl - lose 



-) 



(3) 



(4) 



[0078] As LLR unit 30 is only interested in the relative values of the elements of 
backward vectors Pn, equations (3) and (4) are [preferably] replaced by: 



^+1 = log 



4+1 = log 



+ C 



and 



(3') 



(4') 



[0079] [Preferably, for] the log MAP method, ftmctions f and f-1 are 
implemented partially or entirely using one or more look up tables (LUT). In a 
preferred embodiment of the invention, a LUT is used to approximate the function 
log(l+elx|) and log(l-el'£|). Singularities are [preferably] handled using any of the 
methods described above regarding the linear MAP method. 

[0080] The number LI of symbol groups in each sub-block used in decoding is 
[preferably] chosen as a compromise between the required amount of memory space 
and the decoding speed and accuracy. A conventional number of symbol groups Ll 
used in the prior art is approximately 100. Although substantially any size Ll of the 
sub-block may be used, in some preferred embodiments of the invention, the size Ll 
of the sub-block is substantially larger than that conventionally used in the prior art. 
Using the methods of the present invention, the required size of memory 34 for a 
given size LJ of the sub-block is significantly reduced relative to the prior art. 
Therefore it is possible to increase the size Ll of the sub blocks, achieving faster and 
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better decoding, without substantially increasing the required size of memory 34. In a 
preferred embodiment of the invention, the decoding is performed on entire blocks, 
thus eUrainating the need to accumulate rblodc 52) the additional symbol groups and 
to calculate (blo£ls 56) 3li+1- Thus, the additional time required relative to the prior 

art for performing the first phase fblock 66) of the decoding stage, is partially 
compensated for. 

[0081] It is noted that the size Ll of the sub-blocks used in decoding is not 
necessarily constant and Ll may change during the decoding. 

The implementation of the present invention does not d^end on the number 
of bits included in the state metrics in the forward and backward vectors, and 
substantially any number of bits may be used according to considerations known in 
the art. It is noted, however, that in some preferred embodiments of the present 
invention, the reduction in storage space achieved by the invention is used to increase 
the number of bits included in each state metric and thus increase the accuracy of 
decoder 20. 

[0082] Fig. 6 is a schematic illustration of the backward reference vectors of a 
large block, in accordance with a preferred embodiment of the present invention. In 
some preferred embodiments of the invention, the preparation phase fblock 66i 
comprises a plurality of sub phases, in order to further conserve memory space. 
[Preferably, instead] Instead of calculating and storing all the SM vectors 104 
between the reference vectors in one shot, requiring short term storage 36 to be of the 
size required to store all of vectors 104. backward calculation miit 28 calculates and 
stores the vectors 104 in a pluraUty of sub phases. For example, if the large block 
comprises a thousand symbol groups and tlie number of vectors 104 between the 
reference vectors 102 is one hundred, during the first sub phase backward calmlatmn 
unit 28 may store every tenth vertor 106. In a second sub-phase of preparation phase 
^^IS^ backward calculatinn imit 28 calculates and stores all the vectors 108 
between vectors 106. Thus, for a large block of a thousand symbols, the required 
storage space in memory 34 comprises only room for about thirty vectors: ten vectors 
104 in long term storage 38, and ten vectors 106 and ten vectors 108 in short term 
storage 36. 
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[0083] It is noted that it is not necessary to store entire SM vectors in memory 34. For 
example, in some embodiments the different states of the SM vectors are related in a 
manner which allows calculating the entire vector Pn from half of vector and a 
complementary half of vector Pn+j. In such embodiments, it is possible to store in 

long term storage 38 complementary parts of consecutive vectors rather than entire 
vectors. 

[0084] For example, using as an example, the convolutional code 7.5 described above 
with reference to Fig. lA, rather than storing for each reference vector p^ all the state 
metrics in the vector, backward calculation unit 28 stores the state metrics Pn(01), 
Pn(00), Pn+l(01)> and 3n+l (H)- As can be seen from equations (1), Pn(lO) and 
Pn(ll) can be calculated from Pn+i(01). and 3n+l (H) and therefore storing the 
above state metrics is equivalent to storing vector Pn (ignoring the additional 
calculations). 

[0085] In some preferred embodiments of the present invention, instead of storing all 
the elements of the reference vectors, only those elements of the vectors which are 
considered significant are stored, for some or all of the reference vectors. 
Alternatively or additionally, only some of the elements of the SM vectors stored 
during the decoding stage (in short term storage 36), are stored. The non-signifxcant 
elements are [preferably] elements which have a state metric value much lower than 
other el^ents in the vector such that the chance of that dement to be in a decoding 
trellis is very slight. 

[0086] It is noted that although the above description refers to -forward decoding 
which requires pre-calculation of the backward SM vectors, the principles of the 
present invention apply equivalently to backward decoding in which the forward SM 
vectors are pre-calculated. Furthermore, although decoder 20 was shown as 
including only a single backward calculation unit 28, more than one unit may be used 
in parallel in order to speed up the operation of the decoder. In a preferred 
embodiment of the invention, badcward cqlmilflfmn unit 28 comprises two sub-units 
which operate in parallel during the decoding stage. [Preferably, while] While a first 
sub-unit performs the second phase (71) of a current step, a second sub-unit performs 
the first phase (b^ 66) of a next step. [Preferably, in] M this embodiment, memory 
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34 includes two short term storage spaces 36 which are used by the two sub-units 
interchangeably. 

[0087] It is further noted that the principles of the present invention apply 
substantiaUy to any decoding method which implements the Bail algorithm, and/or 
approximations of the Bahl algorithm, including the linear MAP version, the log MAP 
version, and the approximated log-max version. The principles of the present 
invention may be used with any finite state system, such as communication systems, 
and signal identification systems(e.g., handwriting, voice and image recognition). The 
communication systems may be used in decoding convohitional codes, trellis coded 
modulation (TCM) codes, BCH codes, or any other finite state codes. Alternatively or 
additionally, the principles of the present invention may be applied in performing 
channel equalization and/or in implementing an Inter symbol Interference (ISI) 
equalizer. 

[0088] It will be appreciated that the above described methods may be varied in many 
ways, including, changing the order of steps, and the exact implementation used. It 
should also be appreciated that the above described description of methods and 
apparatus are to be interpreted as including apparatus for carrying out the methods 
and methods of using the apparatus. 

[0089] The present invention has been described using non-Umiting detailed 
descriptions of preferred embodiments thereof that are provided by way of example 
and are not intended to limit the scope of the invention. Variations of embodiments 
described will occur to persons of the art. Furthermore, the terms "comprise." 
"include." "have" and their conjugates, shall mean., when used in the claims, 
"including but not necessarily limited to." The scope of the invention is limited only 
by the following claims: 
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ABSTRACT 

[0090] A method of calculating recursive state metric vectors of a block of symbols 
wherei n , the state metric vectors being supplied to an output calculating unit. The 
method includes [pre-]calculating the state metric vectors from a jSrst side of the 
block, and calculating for a second time at least some of the state metric vectors from 
the first side of the block. 
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